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Abstract. In this note we address the issue of hydrodynamical instabilities in Astrophysical rotating shear flows in the light 
of recent publications focused on the possibility for differential rotation to trigger and sustain turbulence in the absence of a 
magnetic field. We wish to present in a synthetic form the major arguments in favor of this thesis along with a simple schematic 
scenario of the transition to and self-sustenance of such turbulence. We also propose that the turbulent diffusion length scale 
scales as the local Rossby number of the mean flow. A new prescription for the turbulent viscosity is introduced. This viscosity 
reduces to the so-called /J-prescription in the case of velocity profiles with a constant Rossby number, which includes Keplerian 
rotating flows. 



Key words, hydrodynamics - instabilities - turbulence 

' 1. Introduction 

' If it is widely accepted that plane shear flows are subject 
, to non-linear instability, the situation concerning their rotat- 
' ing counterpart is more controversial. Indeed, rotation alone 
, is known to introduce constraints that could result in the sta- 
bilization of the flow. Differentially rotating flows are ubiqui- 
tous in astrophysics. From stellar interiors to accretion disks, 
from galactic rotation to planetary disks, they are met every- 
where. While in many communities it is generally accepted that 
differential rotation could give rise to turbulence, the debate 
, continues in the accretion disk community. Since the dawn of 
■ accretion disk modeling in the early 1970's, and in the ab- 
sence of quantified transport properties supported by an identi- 
fied physical mechanism, turbulence transport has been mod- 
eled with ad hoc prescriptions for the anomalous viscosity 
(the most popular being the so-called "alpha prescription"). 
Differential rotation has been proposed early as the possi- 
ble so urce of this anomalous transport (iShakura & SunvaevI 
Recent numerical simulations fail to show evidence of 
such instability in accretion disks. However, whether or not 
these simulations are adequate or sufficient to give a defini- 
tive answer to this issue rem ains open to question ( Lonaaretti 
l2002t iRichardlEooH l2003h . The difficulty of studying turbu- 
lence and instabilities has been reflected in the lack of a de- 
tailed physical mechanism along with dynamical and transport 
properties to apply in Astrophysical models. Several relevant 
instabilities are currently being investigated in the accretion 
disk community ( Klahr & Bodenheimer (2003), Longaretti 
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(l2002l) . iBalbus & HawlevI (|l99l')), the most studi ed being 
the M agneto-Rotational Instability (MRI) ( Balbus & HawlevI 
Il99ll) . This instability, linear in nature, is an excellent can- 
didate as long as the disk is ionized enough, but is more un- 
likely to be efficient in cold disks such as proto-planetary neb- 
ulae (leading to the int roduction of the "dead zone" layered 
accretion disk model by Gammie (1996)). The issue of differ- 
ential rotation has been recently addressed with rather different 
approaches in several ptiblication such as IChagelishvili et alJ 
('2003), Longaretti' ( 2002), Longaretti' ('20031). iRichard (l2003h 
and Tevzadze et al. ( 2003 ). It is our purpose to try to present 
these different results in a synthetic manner (section |3 and to 
propose in their light a possible basic coherent scenario for the 
development of instabilities and the properties of the resulting 
turbulent state (section|3}- 

2. Flow stability 

In this section we briefly comment on various published results 
concerning the stability of differentially rotating flows. 

2. 1 . The effect of rotation 

If plane shear flows are known to bear non-linear instabilities, 
the picture is quite different where the constraint of rotatio n is 
added. In a widely refeiTed-to analysis, Balb us et alJ d 19961) ar- 
gued that hydrodynamic turbulence would be difficult to sustain 
in Keplerian flows without taking into account hydro-magnetic 
effects. Based on the equations of evolution of the turbulent 
velocity fluctuations energy, they deduced that rotating flows 
could not experience the turbulence due to the energy sink in- 
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troduced by the Coriolis acceleration. Taking a closer look, 
one should note that the Coriolis force embodying the rota- 
tion influence on the fluid motion does not have any effect on 
the total energy equation. Its effect is to redistribute energy be- 
tween the radial and azimuthal components of the motion, van- 
ishing when these equations are summed. This is the natural 
consequence of the fictitious character of the Coriolis force. 
Therefore, the rotation merely stiffens the system, add con- 
straints, but does not introduce any energy sink. The rotation 
and curvature of the flow introduces fundamental diff'erences in 
comparison with plane shear flow. It does not necessary mean 
that the instability mechanisms present in plane flows disap- 
pear, but more likely that the stability and the characteristics 
of the re sulting turbulence are modified (One can refer for ex- 
ample to lLongarettil (l2003l) for a comparative discussion on the 
role of rotation, curvature and viscous dissipation in plane and 
rotating shear flows). 

2.2. Finite amplitude perturbations and transition to 
turbulence 

In contrast with their linear counterpart, non-linear instabilities 
require the presence of fini te (non-infinitesimal) ampU tude per- 
turbatio ns to be tr iggered. Chag ehshvili et all (EooJ) (and se- 
quel Jevzadze.etaL, (i2003)) proposed the mechanism of tran- 
sient growth and bypass mechanism (used for many years in 
the aerodynamics community) to provide the required ampli- 
tudes. Even in a linearly stable flow, the linear operator can 
prov ide transient ampUfication of perturbations. The results 
from'Chageli shvili et alJ ( l2003h show that the initial perturba- 
tion amplitude can be greatly amplified. This amplification can 
introduce into the flow perturbations of ampUtudes sufficient so 
tha t their no n-linear interactions can no longer be neglected. 

iRichardI f2003i) recently proposed a simple model describ- 
ing the necessary conditions for self-sustained turbulence in 
differentially rotating flows, where it can be seen that the en- 
ergy extraction is directly proportional to the shear present in 
the base flow (a classical result for shear flows). It also implies 
that there is a critical Rossby number (Ro - r5,Q/2f2, where 
r is the local radius and Q. is the mean rotation.) above which 
the energy extraction is sufficient to compensate for the stiff- 
ness introduce d in the syste m by the mean rotation. It has to 
be pointed (as IChagelishvili et ak (.20031.) also suggested) that 
the non-linear interactions (also referred to as turbulent diffu- 
sion) do not participate in the energy extraction, but only redis- 
tribute it and counteract the effects of rotation. In this scenario, 
once the critical Rossby number has been reach (Richard 2003) 
and the critical amplitude is present, the flow can then undergo 
a transition from its laminar state to a state where non-linear 
shear turbulence is developed. 



3. Turbulent flow properties 

iRichardI ( l2003h focused on the stability properties of differen- 
tially rotating flows. Using a similar approach, we wish here 
to address the properties of the turbulent state of the bifurcated 
flow, in particular the turbulent transport and scale properties . 



3.1. Rossby number and turbulent scaling 

Rotation has a known tendency to constrain larger scales in 
a two-dimensional state while efficient turbulent diffusion is 
achieved by three-dimensional motions. Typical geophysical 
rotating flows exhibit both two-dimensional (at large scales) 
and three-dimensional (at smaller scales) structures. Within this 
picture, the relevant scale determining the turbulent diffusion 
should be the larger scal e for which the turbu lence can achieve 
to be three-dimensional lBaroud et alJ ( l2003l) performed labo- 
ratory experiments on a rotating annulus, and have shown that 
low Rossby numbers are associated with two-dimensional tur- 
bulence whereas higher Rossby numbers (of order unity in their 
experiments) are associated with three-dimensional turbulent 
structures. The Rossby number for the mean flow can be re- 
written as 



1 r 
Ro= T, 

2 (5,lnQ)-i 



(1) 



where it appears as the ratio of local radius over the char- 
acteristic length scale of the shear. It can also be linked to a 
quantity more often referred to in astrophysics, the epicychc 
frequency w, through the relation 



4Q2 



= (1 +Ro). 



(2) 



The Rossby number of a turbulent structure of character- 
istic length scale A and velocity u ("turbulent Rossby num- 
ber" hereafter), rotating with mean flow velocity, can be ap- 
proximated by Ro^ oc rujlA^Q. : The denominator of the 
Rossby number is twice the rotation experienced by the tur- 
bulent structure, i.e. 2{Q. + u/r) ^ 2Q ; The numerator is the 
derivative with respect to the radius of the turnover time, i.e. 
r5rO + rdriujA) ^ rd,il + ru/A^ ^ ru/A^, where we have used 
the relation 



u/Aoc rd,Sl, 



(3) 



from lRichardl(l2003l) . and the condition A « r. We can 
then write, from Eq.(|3} that 



Roa oc —Ro. 
A 



(4) 



Hence, the ratio of the characteristic turbulent length scale 
over the local radius writes as 



A Ro 
7 R(u' 



(5) 



From Eq. 0, it follows that (for a given Rossby number 
associated with the mean flow) the turbulent Rossby number 
increases when going to smaller scales along the turbulent cas- 
cade. Remembering that high (resp. low) Rossby number are 
associated with three- (resp. two-) dimensional motions, this 
result is coherent with the classic picture of a turbulent spec- 
trum showing two-dimensional structures at large scales and 
three-dimensional ones at smaller scales. 
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3.2. Turbulent diffusion 

From an Astrophysical point of view, the main interest in study- 
ing turbulence might well be to characterize the turbulent trans- 
port of angular momentum or mass, by deriving the relevant 
model for the turbulent viscosity, which is classically written 
as 



V, oc U ■ Ay 



(6) 



where Ay stand for the largest turbulent scale participating 
in the diffusion process and u its associated velocity. From Eqs. 
and (|5jl, it derives that 



i Ro 



PdrQ.. 



(7) 



Two of the results from Richard & Zahn ( 1999), derived 
from experimental data from .Wendt. til933.) and .Tavlor. (.1936) . 



- The ratio Ayjr is constant. From Eq.(|5} it follows that 
(Ro/Ro^J is a constant. 

- The turbulent viscosity scales as v, = y6 • r^drO., where 
y6 is a numerical constant. This result combined with Eq.0, 
implies that /3 oc {RolRox^)^- 



Ex periments from lTavloJ ( Il936l) used bv iRichard & Zahi] 
(11999') to derive the p formulation of turbulent viscosity in ro- 
tating flows, all had the same maximum Ro. The reason is that 
the Rossby number within a laminar Couette-Taylor flow when 
the inner cylinder is at rest is given by Ro - Ri/r, where Rj is 
the radius of the inner cylinder and r is the local radius. It fol- 
lows that the absolute value of the Rossby number is maximum 
at the inner boundary of the flow and decreasing outward. In his 
classic experiments, Taylor modified the aspect ratio of his ap- 
paratus by changing the radius of the outer cylinder, hence did 
not modify the maximum Rossby number value in the flow. It 
follows that in Eq. (Q, Ro is to be treated as a constant for this 
set of experiments. We conclude that both the apparent scal- 
ing of Ay and the constant value of the /3 parameter derived in 
iRichard & Zahn C1999) . can be explained if we postulate that 
there exists a critical (transitional) Ro,i, in the turbulent cas- 
cade above which the turbulence is three-dimensional and un- 
der which it is two-dimensional, and that its value is a constant 
that does not depend on the flow itself. Concluding on the form 
of the turbulent viscosity, in the general case (from Eq.(0) it 
writes 



1 (r25,Q)3 



Ro, 



In the particular case of Keplerian flows, it reduces to 

V, oc ■ r^, 



(8) 



(9) 



where Q.^ = -\/GM/r^ is the Keplerian angular velocity. 
This form is similar to the /3-viscosity prescription, first in- 
troduced by .Lvnden-Bell & Pringle ( ,1974) and more recently 



reintroduced bv IRichard & Zahij (Il999l) (base d on laboratory 
fluid exper i ment data) an d Duschl et al. {20 0^, and apphe d by 
Hure et al.' ("2001"), 'Davis' (i2003i) . iWehrstedt & Gaill (l2003l) and 
Lachaume et aL (.2003.) . 



4. Discussion 

Based on recent developments, we have proposed a simple 
scheme for the transition and self-sustenance of turbulence in 
non-magnetic differentially rotating flows. In this scenario, ini- 
tial velocity fluctuations are amplified by their interaction with 
the mean shear (provided by the linear operator in the equations 
of motion). While this coupling between perturbations and av- 
eraged flow can not by itself destabilize the flow (due to its 
linearly stable nature) it can provide transient growth introduc- 
ing finite amplitude fluctuations. These fluctuations in turn feed 
the non-linear interactions and trigger the instability, under the 
condition that the Rossby number of the mean flow has reached 
a critical value sufficient for the non-linear energy transfer to 
overcome the stiffness introduced by the mean rotation. We 
conjecture that once the flow has become unstable, two- and 
three-dimensional turbulence coexist at different scales. The bi- 
or three-dimensional character of the motions is dictated by the 
value of the Rossby number at a given scale. Finally we have 
introduced a new formulation for the turbulent viscosity of such 
turbulence, which reduces to the /3 prescription in the case of 
constant Rossby number flows, including Keplerian rotation. 

Acknowledgements. The authors wish to thank Dr. Jeffrey Cuzzi and 
Dr. Robert Hogan for their comments and support. D.Richard is 
supported by a Research Associateship from the National Research 
Council / National Academy of Sciences. This research has made use 
of NASA's Astrophysics Data System. 



References 

Balbus, S., Hawley, J., & Stone, J. 1996, ApJ, 467, 76 

Balbus, S. A. & Hawley, J. F 1991, ApJ, 376, 214 

Baroud, C. N., Plapp, B. B., Swinney, H. L., & She, Z. 2003, 

Physics of Fluids, 15,2091 
Chagelishvih, G. D., Zahn, J.-P, Tevzadze, A. G., & 

Lominadze, J. G. 2003, A&A, 402, 401 
Davis, S. S. 2003, ApJ, 592, 1193 

Duschl, W. J., Strittmatter, P A., & Biermann, P L. 2000, 

A&A, 357, 1123 
Gammie, C. F 1996, ApJ, 457, 355 

Hure, J.-M., Richard, D., & Zahn, J.-P 2001, A&A, 367, 1087 
Klahr, H. H. & Bodenheimer, P 2003, ApJ, 582, 869 
Lachaume, R., Malbet, F, & Monin, J.-L. 2003, A&A, 400, 
185 

Longaretti, P 2002, ApJ, 576, 587 
— . 2003, http://www.ArXiv.org, physics/0305052 
Lynden-Bell, D. & Pringle, J. E. 1974, MNRAS, 168, 603 
Richard, D. 2001, These de Doctorat (Universite Paris VII, also 

available at http://tel.ccsd.cnrs.fr) 
— . 2003, A&A, 408, 409 
Richard, D. & Zahn, J.-P 1999, A&A, 347, 734 
Shakura, N. I. & Sunyaev, R. A. 1973, A&A, 24, 337 



4 



Richard & Davis: A note on differentially rotating flows 



Taylor, G. 1936, Proc.Roy.Soc.London A, 157, 546 
Tevzadze, A. G., Chagelishvili, G. D., Zahn, J.-R, & 

Lominadze, J. G. 2003, A&A, 407, 779 
Wehrstedt, M. & Gail, H.-R 2003, A&A, 410, 917 
Wendt, R 1933, Ing.Arch., 4, 577 



